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Abstract. We point out that the nonlinear Schrodinger lattice with a saturable nonlinearity also 
admits staggered periodic as well as localized pulse-like solutions. Further, the same model also 
admits solutions with a short period. We examine the stability of these solutions and find that the 
staggered as well as the short period solutions are stable in most cases. We also show that the 
effective Peierls-Nabarro barrier for the pulse-like soliton solutions is zero. 
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The saturable discrete nonlinear Schrddinger equation is increasingly finding applications in 
various physical situations. Most notably it serves as a model for optical pulse propagation in 
optically modulated photorefractive media [lj, and in this context the pulse dynamics it describes 
have been intensely studied [21 El H]. In addition to its important role for such applications the 
saturable discrete nonlinear Schrddinger equation is also of interest from a purely nonlinear science 
viewpoint [3 El [7]. This interest arises because the saturable discrete nonlinear Schrddinger equation 
has been demonstrated [8] to admit onsite and intersite soliton solutions, which have the same 
energy. This is contrasted by the standard cubic nonlinear Schrddinger lattice where the onsite 
solution always has lower energy than the intersite solution. This phenomenon have often been 
characterized in terms of a so-called Peierls-Nabarro (PN) barrier, which is the energy difference 
between these two distinct solutions. The particular feature of the saturable discrete nonlinear 
Schrddinger equation is thus that it allows the PN barrier to change sign and specifically vanish for 
certain solutions. The vanishing of the PN barrier have been associated with the ability of these 
solutions to translate undisturbed through the lattice, which is impossible in the cubic discrete 
nonlinear Schrddinger equation. Here we derive analytical solutions to the saturable discrete 
nonlinear Schrddinger equation and demonstrate that the localized soliton solutions have a zero 
PN barrier. 

Recently, we obtained |9j two different temporally and spatially periodic solutions to the 
saturable equation [TO] 



i + n\ip n r 

where ip n is a complex valued 'wavefunction' at site n, while v and \i are real parameters. In 
particular, the first solution is 






(2) 
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where the modulus of the elliptic functions m must be chosen such that 

v 2dn(/3,m) 2K(m) 
1 ^"/i" cn2(/3,m)' ^" iV p ' ldj 

and c and 5 are arbitrary constants. We only need to consider c between and | (half the lattice 
spacing). Here Kim) denotes the complete elliptic integral of the first kind [11J . The second solution 
is 

Vff = W- dn( ^ m) cn([n + c]0, m) exp (-»[u;t + 5\) , (4) 

where the modulus m is now determined such that 
_ i/ _ 2cn(/3,m) _ 4ftT(m) 

^ — — — — 2 -, p — — . [0) 



V 



The integer N p denotes the spatial period of the solutions. In the limit N p — > oo (m — > 1), both 
the solutions tpi an d ipn reduce to the same localized solution 

where /3 is now given by 

2-co = - = 2cosh/3 . (7) 
fi 

In Ref. P|, we also developed the stability analysis and examined the linear stability of these 
solutions to show that the solutions are linearly stable in most cases. 

The purpose of this note is to point out that the same model (CE]) also admits the corresponding 
staggered solutions. In particular, using the identities for the Jacobi elliptic functions [12], it is 
easily shown that the model admits the following solutions 

V>f = (-l)^^§^dn([n + c]/3, m) exp (-i[ut + 5]) , (8) 

where the modulus m must be chosen such that 

u o ^_ 2dn(/3,m) ^ ^ _ 2K(m) ^ 
/i cn 2 (/3, m) ' N p 
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€ W = (-l)\l^p^cn([n + c]P,m)exp(-i[u;t + 5}) , (10) 
y fi dn(p, m) 

where the modulus m is now determined such that 

u-2 = -^ = 2cn 2 ^ m ], P = ^M. (11) 



H dn z (/3,m) N } 



p 



In the limit N p — > oo (m — > 1), both the solutions ipi S an d ^n S reduce to the same localized 
staggered solution 

yTi cosh([n + c\p) 

where /3 is now given by 

u-2 = -- = 2cosh/3. (13) 

As an illustration we have plotted the exact solutions of the type IS and IIS in Fig. [TJ Here the 
period N p has to be even. We have shown two periods for type IS and only one for type IIS. 

There are, as expressed by Eqs. (J2]), f TTTT) . and (TlBl . stringent conditions on the parameters 
and z/ for which these exact solutions exist. For example, while the nonstaggered solutions are only 
valid for v > and hence u < 2, the staggered solutions are valid only if v < and hence u > 2. 
In the case IS the limitation is 

(m = 1) < -— < cos 2 f — ) (m = 0) , (14) 

while in the case IIS the limitation is 

( m = 1) < J± < (m = 0) . (15) 

Similarly, the solution ip^ IS ex ists only when — ^ is close to zero (m=l). 

We have also examined the linear stability of these solutions and find that the solutions are 
linearly stable in most cases. A single period (N = N p i, where iV is the lattice size) is always stable 
for both solutions IS and IIS. A type IIS solution with more than one period (N = jN p , where j 
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Figure 1. Illustration of the exact solutions of two types. The parameters are: v = — 1, /i = 0.4, 
cj = 4.5, and c = t = 5 = 0. N p = 10 for and 20 for ipn IS ■ The dashed curves represent the 
solutions given by Eqs. (jHJ and (fpQ|) as if n is a continuous variable. Lines are guides to the eye. 



is an integer larger than one) is also stable, while a type IS solution with more than one period is 
always unstable. Thus, the first example in Fig. [I] is in fact unstable. 

For the solution HIS, expressions for both the power and the Hamiltonian are identical to those 
for the solution III and are given by Eqs. (13) and (14) of Ref. [9]. Hence the PN barrier for the 
solutions III and HIS is the same. We would like to point out here that the calculation of PN barrier 
in I was not quite correct. In particular, since both power P and the Hamiltonian H are constants 
of motion, one must compute the energy difference between the solutions when c = and c = 1/2 
in such a way that the power P is same in both the cases. On using the expressions for P and H as 

given by Eqs. (13) and (14) of I, we find that H for the solution III as well as HIS is given by 

4sinh(/3) 2Bu ( v \ 

H = — + —— + 2 I 1 — — I P . (16) 

/i fl z \ 2/J, J 

Note that H is in fact independent of c, i.e. contrary to our claim in Ref. j9], the PN barrier is in 
fact zero for our solution III (and hence also for HIS). 

Before completing this note, we would like to mention that the model (JT|) also admits a few 
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short period solutions. 
Using the ansatz 

Mt) = <P n e' l{uJt+5) , (17) 
in Eq. ([1]) it is easily checked that the only possible short period solutions are to Eq. ([1]) are: 

(i) Period 1 solution <f) n = (..., a, a, ...) provided 

2 

ua , . 

u ) = -—— i . (18) 
1 + /icr 

(ii) Period 2 solution (p n = (...,a, —a, ...) provided 

2 

ua , . 

w = 4-— — . 19 
1 + yucr 

(iii) Period 3 Solution n = (..., a, 0, —a, ...) provided 

2 

z/a . . 

w = 3 - — ■ — 20 
1 + fia z 

(iv) Period 4 Solutions (fi n = (..., a, 0, —a, 0, ...) and (..., a, a, —a, —a, ...) provided 

2 

z/a , , 

u = 2-—— 5 . 21 

(v) Period 6 Solution n = (..., a, a, 0, —a, —a, 0, ...) provided 

1 + /ia z 

Applying the stability analysis developed in Ref. [9] we have examined the stability of these short 
period solutions and find that for a small nonlinearity < 2fi) they are all stable. The period 
4 solution (..., a, a, — a, — a, ...) is always stable while all the other short period solutions possess 
regions of instabilities at larger nonlinearity. For these low period solutions the stability matrices 
given by Eqs. (20) and (21) of Ref. [S] are simple and it is, for example, easy to see that the lowest 
non-zero eigenvalue, Ai(a, v\ of the stability problem for the period 1 solutions is given by {ji = 1) 

Ai(a,i/) = a 4 + (2-^a 2 + 1. (23) 



Staggered and short period solutions of the Saturable Discrete Nonlinear Schrddinger Equation 7 
Similarly, we have for the period 2 solution 

A 2 (a, v) = a 4 + (2 + a 2 + 1, (24) 

and the period 4 solution 

A 4 (a,z/) = a 4 + (2 - \u\)a 2 + 1. (25) 

These expressions correspond to the relevant curves in Fig. 2. 

It possible to derive similar expressions for the period 3 and period 6 solutions but the 
expressions are more complicated and will be omitted here. Clearly the p period solutions are 
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Figure 2. Regions of stability for the short period solutions to Eq. (TTJ) for /j, = 1. Period 1: thick 
full curve, period 2: dashed-dotted curve, period 3: long-dashed curve, period 4: thin full curve, 
and period 6: short-dashed curve. The instability occurs in the parameter region encompassed by 
the respective curves. 

unstable for the parameter values where A p (a, v) < 0, and we have illustrated these regions in Fig. 
|2j Figure |2] shows the curves in the (a, z/)-plane where X p (a, v) = so that the instability occurs in 
the regions that are encompassed by the respective curves. A symmetry is apparent in this stability 
diagram and it is easy to realize that this symmetry arises from the fact that the transformation 
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(u,(f) n ) — > {—v, ( — l) n 0n) establishes the following connection between the short period solutions: 1 
•H- 2, 3 <H- 6, and 4 <-> 4. 

In conclusion, we have obtained staggered as well as short period solutions of the saturable 
discrete nonlinear Schrddinger equation. We also studied the linear stability and found the solutions 
to be stable in certain parameter ranges. Finally, we found that the Peierls-Nabarro barrier for the 
pulse solutions is zero. Our results are relevant to optical soliton pulse propagation in waveguides 
and photorefractive media [T|. 
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